PROBABILITY CALCULATIONS UNDER THE IAC HYPOTHESIS 
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Abstract. We show how powerful algorithms recently developed for counting lattice points and com- 
puting volumes of convex polyhedra can be used to compute probabilities of a wide variety of events 
of interest in social choice theory. Several illustrative examples are given. 

1. Introduction 

Much research has been undertaken in recent decades with the aim of quantifying the probability of 
occurrence of certain types of election outcomes for a given voting rule under fixed assumptions on the 
distribution of voter preferences. Most prominent among these outcomes of interest are the so-called 
voting paradoxes, which have been shown to be unavoidable, hence the interest in how commonly they 
may occur. The survey [14J discusses these questions and gives a summary of results up to 2002. 

In very many cases, particularly under the IAC hypothesis on voter preferences, the calculations 
involved amount simply to counting integer lattice points inside convex polytopes. In the social 
choice literature, two main methods have been used to carry out such computations. The first, 
dating back several decades, decomposes the polytope into smaller pieces each of which can be treated 
by elementary methods involving simplification of multiple sums. This method works fairly well for 
simple problems but requires considerable ingenuity and perseverance to carry out even for moderately 
complicated ones. More recently, more powerful methods have been introduced in [151 [12] but there 
are several recent instances where even these methods did not suffice to solve natural questions about 
3-candidate elections. 

The purpose of the present paper is to point out that there is an established mathematical theory 
of counting lattice points in convex polytopes (and the closely related issue of computing the volume 
of such a region), which has been partially rediscovered by workers in social choice theory. The area 
has recently been the subject of active research (see [I] for a good summary). Several more efficient 
new algorithms have been devised and implemented in publicly available software. 

We aim to apply these new methods to answer questions in voting theory that have proven beyond 
the reach of previous authors. In addition we corroborate, correct, and unify the derivation of some 
previously published results by using this methodology. We believe that the solution of many hitherto 
difficult problems can now be relegated to a trivial computation. This should open the way for social 
choice theorists to tackle more difficult and realistic problems. We note that Lepelley, Louichi and 
Smaoui |21] have recently, and independently from us, circulated a preprint with a similar goal, which 
covers very similar ground. The fact that two groups of researchers discovered this approach almost 
simultaneously shows that the time has indeed come for these methods to be assimilated by the social 
choice community. 

The basic idea is that many sets of voting situations that are of interest can be characterized by 
linear equations and inequalities. The variables are usually the numbers of voters with each of the 
ml possible preference orders, where m is the number of alternatives. The set of such (in) equalities 
defines a convex polytope in R rf for some d, given by Ax < b for some matrix A (here d < ml and the 
inequality may be strict, because we may first use equality relations to eliminate variables and reduce 
dimension). Each lattice point will correspond to a voting situation in the desired set. The probability 
that a randomly chosen situation has the property under consideration is therefore a straightforward 
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ratio of lattice point counts. Dividing through by n, the total number of voters, yields a convex 
polytope P, independent of n, in For a given number n of voters, the dilation nP describes the 
set of lattice points that we wish to enumerate. 

2. Counting lattice points in convex polytopes 

We give only a brief description here. For more information we recommend [I]. 

The Ehrhart series of the rational polytope P is a rational generating function F(t) = P(t) /Q{t) = 
a n t n whose nth Maclaurin coefficient a n gives the number of lattice points inside the dilation nP. 
The function / : n i— >■ a n is known to be a polynomial of degree d if all the vertices of P are integral; 
otherwise it is a quasipolynomial of some minimal period e. That is, the restriction of / to each fixed 
congruence class modulo e is a polynomial. 

It is known that e is a divisor of m, where m is an integer such that all coordinates of vertices of mP 
are integers. The least such m is the least common multiple of the denominators of the coordinates 
of the vertices of P when each coordinate is written in reduced terms. However there are examples 
where e < m [23 i. A method for determining e was presented in |15j . 

Many questions in voting theory are of most interest in the asymptotic case where n — > oo. For 
small n, issues such as the method of tiebreaking used assume great importance, whereas in the limit 
such issues disappear (the situations in which ties occur correspond in the limit to the boundary of 
P). We focus on limiting results in the present paper. 

The leading coefficient of the quasipolynomial / is the same for all congruence classes: only the 
lower degree terms differ. It is well known that this leading coefficient is precisely the volume of P. 
For many purposes, knowledge of this coefficient is sufficient. The limiting probability under I AC as 
n — > oo is simply the volume of P divided by the volume of X where X is the analogously defined 
polytope that describes all possible voting situations. 

To compute the number of lattice points in nP, if that amount of detail is desired, we may use 
one of several algorithms. An attractive approach pioneered by Barvinok makes heavy use of rational 
generating functions; this is implemented in the software LattE [16] . There are also several algo- 
rithms available for volume computation; see [1] for a survey of algorithms, a hybrid of which has been 
implemented in vinci |26j for floating point computation only. One of these algorithms has been used 
in the Maple package Convex [3j, and uses exact rational arithmetic. 

The software LattE gives the Ehrhart series as standard output. In order to extract the quasipoly- 
nomial formula for f(n) from the Ehrhart series, we may use interpolation. On each congruence class 
modulo e, we must evaluate /(n) at d+1 distinct values of n in this class. Given the explicit expression 
F(t) = P(t)/Q(t) and a computer algebra system, such evaluations are trivially obtained (the a n sat- 
isfy a linear recurrence relation with constant coefficients). The Lagrange interpolation formula then 
yields the desired formula for the particular polynomial that is applicable for the given congruence 
class. 

Another (generally less efficient) method of extraction is to decompose F(t) into partial fractions. 
Note that F(0) = 1 and we can arrange so that Q(t) factors as Y\j(^ ~ a jt) f° r some complex numbers 
otj, possibly not distinct. We then have the partial fraction decomposition 

fW = EEMl -cd)~ k 

a k 

where a runs over the roots of Q and k runs from 1 to the multiplicity of a. 

This shows how the periodicity occurs: the factorization of Q(t) will introduce complex roots of 
unity and the terms corresponding to powers of these will simplify on each congruence class. In fact, 
on extracting the coefficient of t n we obtain 
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and the terms a n simplify on each equivalence class modulo e. Note that e = 1 (that is, f(n) is a 
single polynomial) if and only if Q factors completely over the rationals. 

Note that since we know a priori that the coefficient of t n is polynomially growing, all a with \a\ ^ 1 
can be ignored, since their contribution must cancel (otherwise we would obtain terms exponentially 
growing or decreasing in n). Unfortunately this observation does not help in the present case, because 
the Ehrhart series has a denominator of the form Yl^l — t ai ), so all the cij above are in fact roots of 
unity. 

In summary, the Ehrhart series contains all information required to solve the problem of count- 
ing lattice points in polytopes parametrized by a single parameter n. The hardest step is usually 
determining the minimal period e. 

3. Examples 

In this section we compute, using the recipe above, a few probabilities under IAC that have been 
considered in the recent social choice literature. We emphasize problems where older methods have 
not yielded an answer, but also check results obtained by previous authors using older methods. Some 
of these earlier results appear to be incorrect. The use of a computer algebra system such as Maple 
[22 j is essential for some of the more complicated examples. 

3.1. Manipulability. We first consider the probability under IAC that a voting situation in a 3- 
candidate election is manipulable by some coalition. Counterthreats are not considered — we assume 
that some group of voters with incentive to manipulate will not be opposed by the other, naive, voters. 
See [23 E] for more discussion of these (standard) assumptions. 

For the classical rules plurality and antiplurality, the answer is known: 7/24 and 14/27 respectively 
[181119]. These results were derived by the earliest methods described above and required considerable 
hand computation. However, for the Borda rule, no such result has been derived even using more 
sophisticated methods. A good numerical approximation to the limit has been obtained. In [B] the 
authors used the method of [15] to obtain bounds on the solution but could not carry out the full 
computation. Using their method requires interpolation, hence computing the first 6e coefficients of 
the Ehrhart series, where e is the minimal period of the quasipolynomial. They showed that e > 48, 
and since they computed these coefficients by exhaustive enumeration, it was not possible to carry out 
the computation to the end (the number of voting situations is of order n 5 ). They estimated a value 
of 0.5025 for the limit. 

However with more powerful tools the answers are easily obtained. We let rti, . . . ,n,Q denote the 
number of voters with sincere preference order abc, acb, bac, bca, cab, cba respectively, and let Xi = rii/n. 
Then ^ Xj = 1 and Xj > 0. We use the linear systems derived for general positional rules in [25J. 
As shown in [25] we may assume without loss of generality that a wins, b is second, and c last in 
the election (this assumption will only affect lower order terms in our resulting quasipolynomial, and 
this is inevitable when different tie- breaking assumptions are made). Thus we must multiply our final 
answer by 6 since we are only considering one of the 3! equally likely permutations of the candidates. 

Plurality. We first consider the plurality rule. We define polytopes Pb, P c , Pbc as follows. Consider the 
inequalities 

(3.1) < x\ + X2 — X3 — £4 (a beats b (sincere)) 

(3.2) < £3 + X4 — X5 — xq (b beats c (sincere)) 

(3.3) < — x\ — X2 + X3 + Xi + xq (b beats a (strategic)) 

(3.4) < — x\ — X2 + 2x3 + 2x4 — £5 + 2x2 (b beats c (strategic)). 

The polytope Pb (the region where manipulation in favour of b is possible) is defined by the inequalities 
(|3.1|) - (|3.4f) . the equality ^ Sj = 1, and the condition that all Xj are nonnegative. Polytope P c is 
obtained by applying the permutation b c, which induces the permutation x\ O X2,X3 -H- Xs,X4 O 
xq, and Pbc = Pb D P c is just given by the union of the two sets of inequalities defining Pb and P c . 
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The software LattE readily computes the Ehrhart series of each polytope. They are 

_ 12t 12 + 24i n + AAt 10 + 56t 9 + ffit 8 + 6At 7 + 63t 6 + AAt 5 + 30i 4 + 14i 3 + 6t 2 + 2t + 1 
b ~ (l-t) 2 (l-i 3 ) 4 (l + i) 4 (l + i 2 ) 3 

_ 8t 12 + 16t n + 26i 10 + 34i 9 + 38t 8 + AOt 7 + 41 + 30t 5 + 20t 4 + 10i 3 + At 2 + 2i + 1 

H r = 



H 



(i-t 4 y(i-ty(i-t 2 )(i + t + t 2 y 

At 8 + 5t 6 + At 5 + At 4 + At 3 + 2t 2 + 1 



be 



(l- t y(l-t 4 y(l + t + t 2 y 

The series we require is therefore 
H := Hi + H c — Hie 

__ 16t 12 + 32£ n + 57 1 10 + 68t 9 + 78t 8 + 7At 7 + 73t 6 + 50t 5 + 33£ 4 + 14t 3 + 6t 2 + 2t + l 
~ (l-t 4 ) 3 (l-t) 2 (l-t 2 )(l + t + t 2 ) 4 ' 

Note that in order to factor the denominator of H completely we require both a cube root and fourth 
root of 1, hence a field extension of degree 12. Thus we expect the period of the quasipolynomial 
f(n) := [t n ]H(t) to be 12. We may determine the polynomial formula for / on each congruence class 
in more than one way, as described in section [2] 

First, we try interpolation. Consider the polynomial expression valid for f(n) when n = mod 12. 
This is a polynomial of degree 5 in n. We compute the values /(12j) for j = 0, . . . , 5 and then 
determine the unique interpolating polynomial of degree 5 determined by these points, via, say, the 
Lagrange inversion formula. The built-in commands in Maple find this polynomial immediately: the 
answer is 

f(n) = n 5 H n 4 H n 3 H n 2 H n + 1 (n = mod 12). 

J V ' 17280 108 32 32 120 v ' 

As a check, we substitute n = 96 into this expression — the correct answer, namely [t 9e ]H(t) = 
4176821, is obtained. Analogous formulae can be obtained in the same way for the other congruence 
classes modulo 12. For example, the result for n congruent to 6 modulo 12 is 

T 1 3 15 61 

f(n) = n 5 H n 4 H n 3 H n 2 H n + 5/8 (n = 6 mod 12), 

7 v ; 17280 108 32 32 60 ' v ; ' 

while that for n congruent to 1 is given by 

7 5 1 4 341 , 5 2 917 209 , . 10 , 

tin) = n H n H n H n n (n = 1 mod 12). 

^ ; 17280 108 5184 36 17280 1296 v ' 

Note that since the number of voting situations is given by ( n ^ 5 ) = (n + 1) • • • (n + 5)/120, and we 
have only counted one-sixth of the manipulable situations, the limiting probability of manipulability 
is 720 times the leading coefficient of /, namely 7/24. This agrees with the results obtained in [18J. 
Note that the expressions for finite n do not agree, probably because of different tie-breaking assump- 
tions yielding slightly different sets of manipulable voting situations. We use random tiebreaking as 
described in [25) , with a winner being chosen uniformly at random from the set of those with highest 
score; the alternative used in many papers breaks the symmetry by breaking ties in favour of a fixed 
but arbitrary order on the candidates. It is clear from the discussion at the beginning of the proof of 
[18} Theorem 2] that the latter tiebreaking method is used in that paper. 

As mentioned in section [21 another method would be to compute the full partial fraction decomposi- 
tion of H over the extension field of Q generated by a primitive 12th root of 1. This can be done easily 
by Maple. However the result is somewhat messy and the ensuing computation involving binomial 
coefficients is certainly no easier than using interpolation, so we omit it. 
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Borda. We now consider the Borda rule. We can attempt an analysis similar to the above (the 
polytopes are defined in a similar manner, and all coefficients lie in {0, ±1, ±2, ±3}), but we run into 
serious complexity issues in this case. 

The Ehrhart series F^, F c , F(, c given by LattE are such that when F := Ff, + F c — F^ c is simplified, 
its denominator is a product of cyclotomic polynomials (minimal polynomials for roots of unity) . The 
corresponding roots of unity required are of orders whose least common multiple is 2520. So we are 
still faced with the major task of computing e. It is still an open problem as to whether there exists an 
algorithm to determine e which runs in polynomial time (in the input size) when the dimension is fixed. 
A polynomial time algorithm to determine whether an integer p is equal to e was presented in [27] , but 
has not been implemented in software as far as we are aware. Of course, we do not need to know the 
exact value of e, and we could assume it to be 2520. In order to determine exact formulae for f(n) in all 
cases by interpolation, we would require the first 15120 values of f(n). Trying this in Maple we obtain 
an overflow error. However it would be possible in principle to compute these using the recurrence 
supplied by the rational form of F. We do not proceed further along these lines, but we indicate how 
the computation would go. Writing P(t) = J2k b kt k , Q(t) = Y^k c ^t k \ F(t) = P(t)/Q(t) = J2 n a n t n and 
comparing coefficients, we obtain b n = X^o<fc<n c kf( n — k). This constant coefficient linear recurrence 
allows us to determine sequentially /(0), . . . , f(r) where r = deg P, and for n > r we have the defining 
recurrence ^o<fc<n c fc/( n — k) = 0. In the present case deg P = 75 and deg Q = 82, so the computation 
would be rather involved. 

However, we can certainly determine the leading term of the quasipolynomial /, namely the volume 
of a certain region. It is convenient to eliminate xq throughout, using the sole equality constraint 
Xi = 1. In other words we look at the projection onto the subspace xq = 0. Since we are dividing 
by the volume of the projection of the simplex the exact scale factor is unimportant. This projection 
is defined by the conditions Xj > and Yli=i x i — 1 — we cau these the standard inequalities. 
The volume in M 5 of this simplex is easily computed to be 1/5! = 1/120. Recalling the factor of 6 
mentioned above, we shall therefore multiply the volume answer obtained below by 720 to compute 
the limiting probability. 

The volume required is given by inclusion-exclusion as vol(i2&) + vol(i? c ) — vol(i?j fl R c ) where Rb, R c 
respectively denote the region for which manipulation in favour of b or c is possible. 
The conditions describing the sincere outcome reduce, after elimination of Xg, to 

(3.5) 2x\ + 3^3 + 2x4 — %5 > 1; (a beats b (sincere)) 

(3.6) 2x\ + 3^2 — X4 + 2x5 > 1; (b beats c (sincere)) 

while the conditions describing the outcome after manipulation amount to 

(3.7) 3xi + 4x2 + 3x5 < 2 (b beats a (strategic)) 

(3.8) xi + 2x2 + 2x5 < 1 (b beats c (strategic)). 

Now Rf, is defined by the standard inequalities and those in (|3.5p - f|3.8|> . Also R c is obtained by 
applying the permutation b o c, which induces the permutation x\ «-> X2,X3 •<->■ X5, X4 <-> xq, and Rb c 
is given by the union of the two sets of inequalities defining Rb and R c . 

The package Convex [3] immediately yields the answer when given this input. The respective vol- 
umes of R b , R c , R bc are 371/559872, 881/6531840, 170873/1714608000 and the required limit is precisely 
132953/264600 w 0.5024678760. 

The large denominators in the fractions above give a clue to the difficulty of this problem. Convex 
also computes the vertices of the polytope. The least common multiple of the denominators of the 
coordinates of the vertices is 72 for Rb, 504 for R c , 1260 for Rf, c . Thus the minimum period e is a 
divisor of 2 3 • 3 2 • 5 ■ 7 = 2520, as we already knew from above. 

3.2. Condorcet phenomena. See the two surveys and recent book by Gehrlein [101 111} [T3] for more 
information about previous work on this topic. 
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(A n p) c 


(B n p) c 


B 


(Pnc) 


b (inC) 


0.88148 [9] 


0.62963 [9] 


0.91111(8] 


0.61775 


0.53040 


0.81821 


0.92282 


0.98113 



Table 1. (Joint) limiting Condorcet efficiencies of the standard positional rules under 
IAC, 3 candidates 



In p3] Gehrlein and Lepelley state "A very large number of studies (probably more than 50% 
of the studies that have been devoted to probability calculations in social choice theory) have been 
conducted to develop representations for the probability that Condorcet's Paradox will occur, and for 
the Condorcet efficiency of various rules, with the assumptions of IC and IAC." 

Condorcet's paradox. Condorcet's Paradox occurs in a voting situation when there is no Condorcet 
winner — that is, no one candidate beats all others when only pairwise comparisons are considered. 
This occurrence is independent of the voting rule being used. To compute its likelihood, we compute 
the complementary event. 

Suppose that we have 3 alternatives a, b, c. Let C be the event that a is the Condorcet winner. This 
yields inequalities that boil down to 

(3.9) 2x\ + 2x2 + 2^3 > 1 (a beats b pairwise); 

(3.10) 2x\ + 2x2 + 2x5 > 1 (a beats c pairwise). 

Let Pc be the polytope defined by these and the standard inequalities. Then Convex yields vol(Pc") = 
1/384, so that Condorcet's Paradox occurs with asymptotic probability 1 — 3 • 51/384 = 1/16 for IAC 
with 3 alternatives. This is of course a known result dating back several decades. 

Condorcet efficiency. Similarly we may compute the Condorcet efficiency of a given rule, namely 
the conditional probability that it elects the Condorcet winner given that this winner exists. For a 
given scoring rule defined by weights (1, A, 0), let X\ be the event that a is the winner when this rule 
is used. Clearly Pr(A A ) = 1/3. 

These conditions describing X\ amount to 

(3.11) xi + (1 + A)x 2 + (2A - l)x 3 + (A - l)x 4 + 2Ax 5 > A (a beats b with rule A) 

(3.12) 2xi + (2 - A)x 2 + (1 + A)x 3 + (1 - A)x 4 + Ax 5 > 1 (a beats c with rule A). 

The Condorcet efficiency of rule A is Pr(X\ n C)/ Pr(C) which equals 3 • 51(16/15) vo\(P\ n Pc)- In 
the special cases A = 0, 1/2, 1 of plurality, Borda, antiplurality, respectively, we obtain 119/135, 41/45, 
17/27. These last three results were obtained long ago by Gehrlein. 

We can consider further intersections of such events. For example, Gehrlein has computed limiting 
results under IC for the conditional probability that rule A chooses the Condorcet winner given that 
Borda does, that Borda does given that rule A does, and that both rules choose the Condorcet winner 
given that it exists. The answers to these questions are easily found for IAC using the above methods 
and are listed in Table [TJ These have not previously been published as far as we are aware (numbers 
in brackets in that table represent citations). In Table[T]we let A\C denote the event that antiplurality 
chooses the Condorcet winner given that it exists, B\(P PI C) the probability that Borda chooses the 
Condorcet winner given that plurality does, etc. These can be computed easily using the events C 
and X\ above. For example, the entry B\(P n C) corresponds to the probability of the event that 
Borda and Condorcet agree given that plurality and Condorcet agree. This is simply the volume of 
the polytope P\ji H Pc H Po divided by the volume of Po n Pc (the factor of 3 cancels out because we 
are computing conditional probabilities via P(E\\E2) = P(E\ n E<i) / P(E\)). 

We consider even more intersections of such events in the next section. 

In [2] the value of A for which the positional rule with weights (1,A,0) is most Condorcet efficient 
was determined. We call this "rule M" for brevity. The optimal value of A is an algebraic irrational 
number given as the root of a polynomial of degree 8 and to 5 decimal places equals 0.37228. The 
corresponding value of the Condorcet efficiency is approximately 0.92546, only slightly more than that 
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plurality 


rule M 


Borda 


antiplurality 


0.0278 [17j 


0.00131 





0.0295 P2j 



Table 2. Limiting probabi 



ity of Borda's paradox under IAC, 3 candidates 



for Borda. To use this particular value of A in computations similar to those above, it is probably best 
to switch to software that performs floating point computations in order to compute volumes. One 
such is vinci. We obtain for example that the joint Condorcet efficiency of the optimal rule and the 
Borda rule equals, to 5 decimal places, 0.89183. 

Borda's Paradox. We finish here by discussing Borda's Paradox. Some rules can elect a Condorcet 
loser, namely a candidate that is beaten by every other when pairwise comparisons are made. The 
probability of this event for plurality and antiplurality has been studied under IAC in [T7], and it 
has long been known to be zero for Borda. The methods in this section can be applied directly, 
since we need only replace the Condorcet winner conditions by the same ones with the direction 
of the inequality reversed. This shows that Borda's Paradox occurs for plurality with probability 
1/36, agreeing with [17J. The corresponding results for Borda and antiplurality are and 17/576, 
corroborating the previous results. The probability that the most Condorcet efficient rule above elects 
the Condorcet loser is, as one might expect, very small. The results are shown in Table [2j 

3.3. When do all common rules elect the same winner? For three-alternative elections, all 
positional voting rules elect the same winner in a given situation if and only if both plurality and 
antiplurality elect the same winner, since the vector of scores is a convex combination of those for the 
two extreme rules. The probability of this event has been investigated under IC but not under IAC 
as far as we are aware. In |24j Merlin, Tataru and Valognes also investigated the probability under 
IC that all positional rules and all Condorcet efficient rules yield the same winner (in this case, all 
scoring runoff rules also yield this same winner) . 

We again suppose that a is the winner. We want to compute the probability of the event P n 
A as described in the previous section. The relevant polytope has 18 vertices and m = 12. Its 
volume is 113/77760 and so the limiting probability that all positional rules yield the same winner 
for 3 alternatives under IAC is 113/216 (this confirms a result in |12j). We could also investigate 
the relationship between, say, plurality and Borda. They agree with probability 89/108, whereas 
antiplurality and Borda agree with probability 1039/1512. 

The probability that all Condorcet rules and all positional rules elect the same winner given that 
the Condorcet winner exists is obtained easily via computation of Pr(PnCn^4) as above. The answer 
is 3437/6480. The polytope involved has 29 vertices and m = 12. 

We must also consider the case when no Condorcet winner exists. There are two cases corresponding 
to the two cycles a, b, c, a and a, c, b, a. In the first case, [21] shows that the rules all agree if and only if 
all positional rules give the ranking a, b, c, and this occurs if and only if both plurality and antiplurality 
give that ordering. The computation is straightforward as above and the probability of this event is 
only 5/10368. The contribution from the cyclic case is therefore 32 times this, or, 5/324, and the final 
result for the probability that all rules agree is 10631/20736. 

We can also consider the probability that two rules agree in their whole ranking, not just in the 
choice of winner. This is easily computed similarly to above: plurality and antiplurality agree on their 
whole ranking with probability 8/27, while Borda and plurality agree with probability 61/108. Borda 
and antiplurality also agree with probability 61/108, which is clear by symmetry in any case. 

3.4. Abstention and Participation Paradoxes. In [20] Lepelley and Merlin discuss various ways 
in which voters can attempt to manipulate an election by abstaining from voting. All scoring runoff 
rules and Condorcet rules suffer from this problem. Although abstaining turns out to be a dominated 
strategy for scoring runoff rules, it is still of interest to compute the probability that a situation may be 
manipulated in this way. Lepelley and Merlin carry this out under IC and IAC for scoring runoff rules 
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Rules 


Elect same winner 


Agree whole ranking 


Antiplurality and Borda 


0.68717 


0.56481 


Antiplurality and plurality (hence all scoring rules) 


0.52315 [12] 


0.29630 


Plurality and Borda 


0.82407 


0.56481 


All common rules 


0.51268 




Table 3. Limiting probability of agreement of various rules under '. 


AC, 3 candidates 



Underlying rule 


PPP 


NPP 


PAP 


NAP 


plurality [2Ti] 
Borda 
antiplurality [20j 




0.01389 
0.03822 


0.07292 
0.02083 





0.01042 
0.04253 


0.04080 
0.01389 




Table 4. Limiting probability of participation paradoxes for scoring runoff rules under 
I AC, 3 candidates 



based on plurality, antiplurality and Borda. For the latter (the Nanson rule) the limiting probability 
was not computed exactly (Table 5 of [20] refers to results of Monte Carlo simulation). We compute 
some exact values here. 

We use the linear system given in [20]. Suppose that c is eliminated first and a then beats b in 
the runoff. The Positive Participation Paradox occurs when voters ranking a first are added to the 
electorate, and yet a then loses. This cannot happen when plurality is used at the first stage, but for 
other rules it can happen that b now loses the first stage, and a subsequently loses the runoff against 
c. Note that only voters with preference order acb can cause this to occur, and it can only occur when 
c originally beats a pairwise. 

The system describing this set of voting situations contains the inequalities stating that a beats b 
and b beats c using the given scoring rule, and also that a beats b pairwise. In addition we have another 
constraint as described in [20] (note that rig in the first equation on p. 58 of that paper should be — n%). 
Carrying out the (by now routine) computation we obtain 1/72 which confirms the simulation result 
0.14 referred to above. Note that the polytope involved has only 6 vertices and 6 facets but m = 18; 
if e = 18 (which we have not checked), it would be difficult to compute the Ehrhart polynomial using 
the old methods, which probably explains why only simulation results were obtained for the Nanson 
rule in the paper cited above. 

Similarly we may compute the result for each of several other participation paradoxes. The results 
for the negative participation, positive abstention and negative abstention paradoxes (see [20] for 
definitions and characterizations of the polytopes) are respectively 1/48,1/96,1/72 confirming the 
earlier simulation results 0.020,0.010,0.14. 

We can also perform the analogous computations for plurality and antiplurality runoff — the results 
confirm those in [20] and are shown in Tabled! 

3.5. The referendum paradox. This gives an example where the variables describing our polytopes 
are slightly different. 

In [7] the referendum or Compound Majority Paradox is studied. In the simplest case there are N 
equal sized districts each having n voters. There are two candidates a and b and voters in each district 
use majority rule to decide which candidate wins each district. The candidate winning a majority 
of districts is the winner of the election; the paradox occurs when this candidate would have lost if 
simple majority had been used in the union of all districts. 

Among other things, the authors of [7\ derive the probability of occurrence for N = 3,4,5 under 
IAC using the older methods and state that they are not able to extend it to iV > 6. Using the 
methods of the present paper it is easy to perform the computations for at least a few more values 
of N. Let rii denote the number of voters voting for a in district i. The relevant set turns out to be 
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number of districts 


3 


4 


5 


6 


7 


9 


probability 


0.125 [7J 


0.02083 [7] 


0.15885 [7\ 


0.04063 


0.20419 


0.26954 



Table 5. Limiting probability of referendum paradox under IAC, 3 candidates 



described (ignoring ties for simplicity) by the union of polytopes of the form 

Tit > N/2 for 1 < i < k (a wins k districts) 
< rii < N/2 for k + 1 < i < N (b wins N-k districts) 

rij < Nn/2 (b wins overall) 

i 

for \N /2\ +1 < k < N — 1. The polytope P n corresponding to k must have volume multiplied by 2(^f) 
to account for the symmetries of the problem. Note that there are (n + 1)^ situations to consider and 
so the leading term in n of the Ehrhart polynomial gives the probability required. 

Doing the analogous computation for N = 3 and N = 4 we obtain (as n — > oo, in other words 
computing the volume of P n /n) results agreeing with |7J. However already for N = 5 we obtain 
61/384 as opposed to their result 55/384. For N = 7 we have 9409/46080. The most complicated 
corresponding polytope in the last case has 36 vertices and 11 facets, whereas for N = 9 it has 91 
vertices and 14 facets. We did not attempt to find the maximum value of N for which our software 
could obtain an answer; the answer was given essentially instantaneously for N = 9. The conjecture 
in [7] that the probability tends to a limit of around 0.165 as N (odd) goes to infinity seems unlikely 
in the light of these results. 

4. Summary and discussion of future work 

We have shown that a wide variety of natural probabilistic questions for 3-alternative elections under 
IAC can be answered by applying standard algorithms for counting lattice points in, and computing 
volumes of, convex polytopes. For 4 or more alternatives the computations are conceptually the 
same but necessarily more complicated. However, the scope for extending results in the 3-candidate 
case to 4 or more candidates is obviously higher than for the older methods, which now appear to be 
completely superseded. One important point to notice is that many algorithms for volume computation 
have running times that are very sensitive to the number of defining hyperplanes and the number of 
vertices. Thus finding the most efficient description of the input system is important. It is certainly 
clear that further progress in this area will require researchers in social choice theory to understand 
in some detail how the fastest algorithms for lattice point counting and volume computation actually 
work. This may even lead to proofs for larger (or general) numbers of candidates when the polytopes 
concerned have a particularly nice structure. 

Many questions naturally arise from our work here. One obvious line of attack is to try to find the 
optimal parameter for 3-alternative scoring rules that minimizes the probability of a certain undesirable 
behaviour occurring. The present authors are already engaged in carrying this out for the case of 
(naive, coalitional) manipulability. Numerical results obtained in |25j show that the answer may well 
be plurality, but this has never been proved. An attack on this problem along the lines of the approach 
in the present paper would require computation of volumes of a polytope whose defining constraints 
depend linearly in a parameter A, and this requires considerable work as shown in [2]. Understanding 
of how to carry out such a computation would help in understanding the variation between positional 
rules. For example, the probability of electing a Condorcet loser is of order 0.03 for both plurality and 
antiplurality, but an order of magnitude smaller near the Borda rule, and as a function of A is very 
flat there. Quantifying this type of variation analytically may show, for example, that it is not worth 
the trouble of replacing Borda by the Condorcet-optimal positional rule. 

Another direction is to consider other probability models. For simplicity here we have not considered 
some common assumptions such as single-peaked preferences and the Maximal Culture Condition. 
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Many computations in these cases reduce to ones identical in spirit to those we have undertaken here. 
More general Polya-Eggenberger distributions would lead to the more difficult issue of integrals of 
nonconstant probability densities over polytopes, but some results may be forthcoming there. 
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